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We solved with a numerical procedure the HTL improved ladder DS equation for the 
retarded fermion self-energy function Er to study the spontaneous generation of fermion 
mass in thermal QCD/QED, and studied the gauge-dependence of the solution within a 
general covariant gauge where the gauge parameter ^ is any constant number. 

With the numerical solutions thus obtained, we found the foUowings; i) The fermion 
wave function renormalization function A{P) always deviates largely from unity even at the 
momentum where the mass is defined, thus the corresponding solutions explicitly contradict 
with the Ward-Takahashi identity, ii) As a result, the obtained solutions strongly depend on 
the choice of gauge parameters: the critical temperatures and the critical coupling constants 
significantly change gauge by gauge. In all gauges we studied in the present analysis, we 
could not find any solution, having a possibility to be consistent with the Ward-Takahashi 
identity. Thus we are forced to investigate the procedure to find a gauge which enables us 
to get a solution being consistent with the Ward-Takahashi identity, otherwise we can not 
obtain any physically sensible conclusions through the analysis of the point- vertex ladder DS 
equation no matter how the gauge propagator gets improved. 

§1. Introduction and summary 

Recent studies on gauge field theory have been reveahng rich aspects of phase 
structure of the matter according to the variation of number density and/or temper- 
ature. However, it is supposed to be hard to obtain more detailed understanding of 
the mechanism of phase transition by means of such theoretical studies, since most of 
them were carried out on the basis of perturbative calculation or of numerical lattice 
simulation. This motivated us to survey the problem of phase transition using the 
Dyson-Schwinger (DS) equation, firstly because the DS equation is derived exactly 
in the field theory and is the fundamental equation to investigate nonperturbative 
phenomena within the framework of field theory, and secondly because we can ob- 
tain successively improved solutions by the successive refinement of the analytical 
approximation to its integration kernel, thus revealing the essential contribution that 
controls the phase transition depending on the temperature/density. 

We started our analysis with the DS equation for the retarded fermion mass func- 
tion to study the spontaneous generation of fermion mass in thermal QED/QCD.^-* 
In the analysis we used an improved ladder interaction kernel obtained analytically 
through the Hard-Thermal-Loop (HTL) resummation procedure, in which the ladder 
kernel is improved by use of the HTL resummed form of the gauge boson propaga- 
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tor. We realized that the results obtained^-' are significantly different from those 
obtained in the preceding analyses with the simple ladder DS equation. '^^ It is worth 
noticing that, in all preceding DS equation analyses including ours, the bare (point) 
vertex ladder approximation for the integration kernel has been used, thus no result 
considering any vertex correction. 

We here summarize the essential points of our analysis that may give the result 
significantly different from those of the preceding works. (For details of our results, 
see Ref. 2).) 

i) The DS equation for the retarded fermion mass function Ufj is derived correctly 
without any specific assumption for its form. 

ii) Nontrivial imaginary parts of the invariant functions A, B and C (see, the 
definition of Eq. (jl-ip ) are taken into account. 

iii) We adopted the improved ladder integration kernel analytically obtained in the 
HTL approximation. 

Despite the improvements in taking into account correctly a) the (unstable) 
thermal quasiparticle nature of the fermion in the heat bath and b) the dominant 
effect of thermal fluctuations through the HTL resummation, it was suggested that 
the obtained result showed a serious problem, i.e., dependence on the choice of the 
gauge used in the analysis. We can recognize it to see that the Ward-Takahashi 
identity Zi = Z2 does not hold, where Zi and Z2 denote the vertex and the wave- 
function renormalization constants, respectively. 

The retarded fermion mass function can be parameterized with the three invari- 
ant functions A, B and C as 

Sr{P) = (1 - A{P))pa' - B{P)j^ + C{P) (M) 

The inverse of A{P) at the momentum where the fermion mass is calculated is 
nothing but the wave function renormalization constant Z2. As noted above, the 
vertex renormalization constant Zi is exactly unity, Zi = 1, in our analysis because 
of the bare (point) vertex ladder approximation for the integration kernel, which is 
also the case with other works carried out before. It follows that the Ward-Takahashi 
identity, the statement of gauge invariance, requires Z2 = 1, namely, A{P) = 1 at 
least at the momentum where the fermion mass is calculated. If ^4 7^ 1 in the 
obtained results, it means that the results do not satisfy gauge invariance and hence 
there is little physical meaning in the results obtained. 

For the vacuum (T = 0) case, fortunately, A{P) = 1 is verified to hold in the 
analysis in the Landau gauge. Therefore it is supposed that the analysis of the DS 
equation in the Landau gauge, apart from the reliability of the ladder approximation, 
has some physical significance in the quantitative as well as qualitative sense. ^^'^^ 
At finite temperatures, however, there is no guarantee that A = 1 holds even at 
the momentum where the fermion mass is calculated. In fact our analysis in the 
Landau gauge shows that A largely deviates from 1 and becomes even complex 
number.^) Namely the analysis of DS equation in the ladder approximation at finite 
temperature/density QED/QCD is obviously inconsistent with gauge invariance. 

This gauge-dependence problem must be tackled seriously in order to draw a 
definite conclusion from our analysis of the HTL improved ladder DS equation. 
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which indicated the importance to correctly take the dominant effect of thermal 
fluctuations into the integration kernel through the HTL resummation. 

To see the problem more clearly we must clarify how sensitive the obtained 
results are to the choice of gauges. It will help us to study whether there exists 
a solution of the DS eqution in the ladder approximation that satisfies the Ward- 
Takahashi identity Zi = Z2 at finite temperature/density, and also to investigate, if 
such a solution exists, what is the real difference of it from the "gauge-dependent" 
solutions. 

To answer the questions, in this paper we will solve the HTL improved ladder DS 
equation in a general covariant gauge and study the dependence of the solutions on 
the various choices of the gauge parameter, then investigate the possibility to obtain a 
solution consistent with the Ward-Takahashi identity. We also make an improvement 
in estimating the numerical integration over the singular part of integration kernel 
in the present analysis. 

We here present the results of our analysis; We find, within gauges where the 
gauge parameter ^ is constant numbers, i) that the fermion wave function renor- 
malization function A{P) always deviates largely from unity even at the momentum 
where the mass is defined, thus that the corresponding solutions explicitly contradict 
with the Ward-Takahashi identity, and ii) that, as a result, the obtained solutions 
strongly depend on the choice of gauge parameters: the critical temperatures and 
the critical coupling constants significantly change gauge by gauge. In all gauges we 
study in the present analysis, we can not find any solution that has a possibility to be 
consistent with the Ward-Takahashi identity. Thus we are forced to investigate the 
procedure to find a gauge which enables us to get a solution being consistent with 
the Ward-Takahashi identity, otherwise we can not obtain any physically sensible 
conclusions through the analysis of the point-vertex ladder DS eqution no matter 
how the gauge propagator gets improved. 

This paper is organized as follows; In the next ^we present the HTL resummed 
Dyson-Schwinger equation for the retarded fermion self-energy function and 
explain the improved ladder approximation adopted in the present analysis. Section 
3 is devoted to presenting the solutions of the HTL improved ladder DS equation, 
revealing the large gauge-dependence between solutions. Conclusions of the present 
paper and the related discussion are given in the last 21 I^i the Appendix we give 
some technical details in the numerical analysis in solving the DS equations. 

§2. The HTL resummed improved ladder Dyson-Schwinger equation 

In this section we present the HTL resummed DS equation for the retarded 
fermion self-energy function We also give an explication about the improved 
ladder approximation we make use of to the HTL resummed gauge boson propagator. 

We then present the HTL resummed improved ladder Dyson-Schwinger equation 
for the independent invariant scalar functions A, B and C. We also calculate the 
effective potential for the retarded fermion propagator Sr in order to find the "true 
solution" when we get several "solutions" of the DS equation. 
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2.1. The HTL resummed Dyson-Schwinger equation for the retarded fermion self- 
energy function 

In the real time closed time-path formalism, we obtain, in the massless thermal 
QED/QCD in the HTL approximation, the DS equation for retarded fermion self- 
energy function Uj^: 



-iURiP) = -iEnAi-P,P) = -y /■ 



(2vr)4 

X [*r^AAi-P^ K,P- K)Sra{K) *r-i^^A-K, P,K-P) *Grr,^,{P - K) 

+ *r^Aj,{-P, K,p- k)Srr{k) *riAR{-K, P,K-p) *GraMp - k)] , 

(2-1) 

Here *G^'^ is the HTL resummed gauge boson propagator, whose retarded (R = RA) 
and correlation (C = RR) components are given by^^ 

1 1 £ 



*n^{K) -K"^ - ieko *n^{K) - K"^ - ieko + ieko 

(2-2) 

*G^c''{K)^*G''^'^{-K,K) = {l + 2nB{ko))[*G%'{K)-*G>'/{K)], (2-3) 

^^Biko) = TTTT^^x 7, (2-4) 

exp(A;o/r) - 1 

with *77^ and *77f being the HTL contributions to the transverse and longitudinal 
modes of the retarded gauge boson self-energy, respectively.*^) The parameter £ is 
the gauge-fixing parameter (£ = in the Landau gauge). In the above, A^'^ , B^" 
and D^'^ are the projection tensors given by"^) 



A^^- = g^^- -B^^- -D^^, B^^ = D^, = ^^, (2-5) 



where K = (k, k^k), k = \/k^ and k = k/A; denotes the unit three vector along k. 

We use S{—P, P) = S{P) to denote the full fermion propagator, whose retarded 
(R = RA) and correlation (C = RR) components are given by 

Sr{P)^Sra{-P,P) = -x—-^ —, (2-6) 

f+ ie-/o - 2Jr 

Sc{P) = Srr{-P,P) = {l-2nF{po))[SR{P)-SA{P)], (2-7) 

npiPo) = . ]rr,^v (2-^) 

exp{po/T) + 1 

with the retarded fermion self-energy function Er decomposed as Eq. (jl-ip in terms 
of the independent invariant scalar functions A(P), B{P) and C{P). 

Finally, the HTL resummed 3-point fermion- gauge boson vertex functions, *r^^, 
are given by 

^RAA = ^AAR = 7^, Otherwise 0. 
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where SF^^^ denotes the HTL resummed contribution to the vertex function.^) 

As mentioned in the introduction, at zero temperature the fermion wave function 
renormahzation constant A{P) is equal to unity in the Landau gauge = 0) even 
in the ladder (point-vertex) DS equation, while at finite temperature it is not. The 
quantity C{P)/A{P) = M{P) plays the role of the mass function, in which we are 
interested, that vanishes in the chiral symmetric phase. 

2.2. The HTL resummed DS equations for the invariant functions A, B and C 

In the present analysis, we solve the DS equation for the retarded fermion self- 
energy function with the HTL resummed gauge boson propagator, by adopting 
further the following two approximations, i) the point-vertex approximation, and ii) 
the modified instantaneous exchange approximation, on which we give brief expla- 
nations below. 

i) Point-vertex approximation 

As for the vertex function *r^^ we adopt the point- vertex approximation, namely 
we disregard SPj^^^ in Eq. (|2-9p . Thus we investigate the ladder (point- vertex) DS 
equation with the HTL resummed gauge boson propagator. 

There are two reasons; Firstly, without the point-vertex approximation the nu- 
merical calculation we should carry out becomes so complicated that we can not 
manage with the power of the computer we use, because the HTL resummed con- 
tribution to the vertex function, drj^^^, is the non-local interaction term, and also 
because it behaves singular in numerical calculations. Secondly, in the DS equation 
with the HTL resummed vertex function, it is difficult to resolve the problem of dou- 
ble counting of diagrams,^''-* especially at the level of numerical analyses. Being free 
from this problem is the main reason we make use of the point- vertex approximation. 

ii) Modified Instantaneous Exchange (MIE) approximation 

We make use of the modified instantaneous exchange (IE) approximation (i.e., 
set the energy component of the gauge boson to be zero) to the gauge boson prop- 
agator, which consists of taking the IE limit in the HTL resummed longitudinal 
(electric) gauge boson propagator, *G'^ , that is proportional to B>^'^ , while keeping 
the exact HTL resummed form for the transverse (magnetic) gauge boson propa- 
gator, *Gjf , that is proportional to A^^" , and also for the massless gauge term in 
proportion to D'^'^ . The reason why we do not take the IE limit to the transverse 
mode is that the IE approximation reduces the transverse mode to the pure massless 
propagation, thus makes the important thermal effect, i.e., the dynamical screening 
of transverse propagation disappears. 

With the above two approximations, we obtain the HTL resummed improved 
ladder DS equations for the invariant scalar functions A^ B and C: 



p\l-A{P)] = e^ j 



(2^ 



{1 + 2nB{po - ko)}Im[ *G'l^{P - K)] 



X 



{K^Pp + KpP„ - poiK^Qpo + Kpg^o) - ko{PagpO + Ppgao) + pkzg. 



crp 



-2pofco5ao5po} 
-2po5f705po} 
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A{K) 



[ko + B{K) + ie]2 - AiKfk"^ - C{KY 
ko + B{K) 



+ {1 - 2nF{ko)} X 



[ko + B{K) + ie]2 - A{K)-^k^ - C{K)^. 
'Cj^iP - K)Im[{K„Pp + KpP^ - poiK^gpo + i^p5ao) - HP^gpO + PpQao) 

A{K) 



+pkzg„p + 2po^o5'aOS'po} 



+{PagpO + PpQuO - 2po5(705po} 



[ko + B{K) + ie]2 - A{Kfk^ - C{Kf 
ko + B{K) 



[ko + B{K) + ief - A{Kfk'^ - C{Kf 



B{P) = e' I 



(2-10) 



d-^K 



(2vr)4 

-{25po25f70 - 5f7p} 



{1 + 2nB(iJo - A;o)}/m[ *Gg^(P - K)] x 



fco + B{K) 



*GP^{P-K)I' 



m 



-2kogaogpo} + 



[ko + B{K) + ie]2 - A(i^)2fc2 - C(K)2 

MK) 

[ko + + ie]2 - A(K)2A:2 - C{KY 

ko + B{K) 



+ {l-2nF{ko)] X 

{K^QpO + -f^pfifcrO 



[fco + B{K) + ie]2 - A(i^)2A;2 - C{KY 



{2gpo2gcro - gap} 



(2-11) 



C{P) 



I ^9ap{l + 2nB{po - ko)}Im[ *G'^{P - K)] 
C{K) 



[ko + B{K) + ief - A{K)^k^ - C{Kf 
GP^{P-K)I-^ ^^^^ 



+ {1 - 2nF{ko)} X 



m 



L [ko + B{K) + ie]2 - ^(i^)2A:2 - C(K)2 



(2-12) 



The HTL resummed DS equations are the coupled integral equations for the six 
unknown functions because the invariants A, B and C have both real and imaginary 
parts. Therefore, the DS equations, Eqs. ()2-10p - (l2-12p . are still quite tough to be 
solved even if we adopt the above two approximations. 

2.3. The effective potential V[Sr\ for the retarded full fermion propagator Sr 

The above DS equations, Eqs. (|2-10p - (l2-12p . may have several solutions, and we 
choose the "true" solution by evaluating the effective potential T^[5'i{] for the fermion 
propagator function Sr, then finding the lowest energy solution. The effective po- 
tential is expressed as^^^ 

VISrHTt IfSs] + iTrlii [iS^'] 
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+1^Sc{K)^uSr{P)D^;^{P -K)+ j^Sr{K)j,Sc{P)D'X'{P - K)] , 

(2-13) 

where the 1st and the 2nd terms correspond to the one-loop contribution, while the 
3rd term corresponds to the two-loop contribution. 

§3. Solution of the HTL improved ladder DS equation and its 

gauge-dependence 

In this section we solve the HTL improved ladder DS equation derived in the 
previous section numericallj0 by an iterative method. To start with we choose ap- 
propriate initial trial functions for A{P), B{P) and C{P) with the guess for them to 
have non-trivial imaginary parts. By calculating the right-hand sides of Eqs. (j2-10p - 
(j2-12p we get "solutions" A{P), B{P) and C{P), which are supposed to be better 
approximations to the real solutions. Replacing the initial trial functions for A(P), 
B{P) and C(P) by the "solutions" thus got, we obtain further better approximations 
to the solutions. This procedure is iterated till we obtain converged solutions. Since 
the solutions thus determined may depend on the initial choice of the trial functions, 
we need to try various initial trial functions to find the real solutions. Some tech- 
nical details in the numerical analysis in solving the DS equations are given in the 
Appendix. 

If more than two converged solutions are obtained, among the solutions we adopt 
such a solution to be the true one that has the lowest value of the effective potential 
VISr], Eq. (EH. 

To investigate how the solution thus obtained depends on the choice of gauges, 
we solve the HTL improved ladder DS equation by choosing various gauges within 
the covariant gauge. In the present paper, we carry out the analysis by choosing 
gauges mainly in the neighborhood of the Landau gauge (.^ = 0), and show explicitly 
the obtained solutions suffer from the strong gauge-dependence. 

Gauges we choose are those with five different values of the gauge-parameter ^, 
including the Landau gauge: 

i = 0(Landau gauge), ^ = ±0.05 and ^ = ±0.025. 

3.1. Result in the Landau gauge = 

Firstly we give the result in the Landau gauge, which is the gauge studied in 
most of the preceding analyses. Fig. 1 shows the mass Re[M[P)] = Re[C {P) / A{P)] 
as a function of temperature T at five different values of the coupling constant a, 
a = 3.5, 3.7, 4.0, 4.5 and 5.0. In this figure we also give the critical exponent v at 
each coupling constant a defined by 

i?e[M] = C7„ (T, - ^)^ r<r„ (3-l) 

which controls how the mass i?e[M] vanishes near the critical temperature Tc- We see 
that the temperature dependence of the mass i?e[M] near the critical temperature 

*-* As is evident from the definition, Eq. the invariants B and C have a dimension of mass. 

Thus, in this section, they are measured in units of the cutoff scale A in all the results shown. 
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Fig. 1. The T-dependence of the mass Re[M] = Re[C/A] at po = 0, p = O.lyl for various fixed 
values of the coupling constant a in the Landau gauge — 0). The best-fit curves at each 
coupling constant, with the critical exponents given in each data, are also shown. 
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Fig. 2. Comparison of the wave function renormalization constant i?e[yl] and |Jm[A]| at the cou- 
pling constant a — 4.0 and 5.0 evaluated at po = 0, p = 0.1 A in the Landau gauge = 0). 



Tc can be well described by the functional form Eq. (j3-ip with little depnedence of 
the critical exponent on the couplng constant. The average value v ~ 0.527 can 
reproduce well the result with appropriate Ca- Thus we may say that the phase 
transition takes place through the second order transitionEl 

*-* To be exact, in order to conclude that the phase transition is of the second order, we must 
study at each coupling constant whether there are no stable states other than the one shown in the 
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Fig. 3. Gauge-parameter-dependence of the fermion mass Re[M] = Re[C/A] at the coupling con- 
stant a = 4.0 evaluated at po = 0, p = O.lyl. The best-fit curves for each value of the gauge- 
parameter, with the critical exponents u given in each data, are also shown. 



However, this result can not be justified in having a physical significance without 
further consideration. As shown in Fig. 2, the invariant function A{P) at the momen- 
tum where the fermion mass is defined deviates largely from 1, even its imaginary 
part being sizable: i2e[A] > 1.4 and |/m[^]| ~ 0.1 ~ 0.2, indicating Z2 signifi- 
cantly smaller than 1, not even the real number. This fact implies that the result 
obtained in the Landau gauge apparently contradicts with the Ward-Takahashi iden- 
tity Zi = Z2, therefore we should not give a serious physical meaning to the result 
in the Landau gauge. 

3.2. Gauge- dependence of the solution 

Now we compare the results with the five different values of ^ in order to see 
how large the solution depends on the choice of gauge. We illustrate it by showing 
the results obtained in the case of the coupling constant a = 4.0. Fig. 3 shows how 
the mass M{P), as a function of temperature T, depends on the choice of gauge- 
parameters. 

In this figure we give also the critical exponent at each value of the gauge- 
parameter. In any gauge the temperature-dependence of the mass Re[M] near the 
critical temperature can be well described by the functional form Eq. (j3-ip with little 
depnedence of the critical exponent on the couplng constant. The average value 
v ~ 0.527 can reproduce well the result with appropriate Cq,. This fact indicates 
that the phase transition takes place through the second order transition. 

However, as can be easily seen, the critical temperature strongly depends on the 
choice of gauge-parameter, demonstrating that it is essential to choose an appropriate 



figure. 
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Fig. 4. Comparison of the wave function renormalization constant |^| at the coupling constant 
a = 4.0 evaluated at po = 0, p = O.lyl. 
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Fig. 5. Comparison of the wave function renormalization constant iie[yl] and |/m[j4]| at the cou- 
pling constant a = 4.0 evaluated at po = 0, p = 0.1/1. 



gauge such that the Ward-Takahashi identity holds in solving the improved ladder 
DS equation in order to obtain reliable results which have predictive power. 

For further consideration of the optimal gauge, we show the gauge dependence of 
the wave-function renormalization function A{P) in Figs. 4 and 5. Fig. 4 shows the 
absolute value of A[P) as a function of temperature for the five gauge parameters, 
and Fig. 5 shows the real and imaginary parts of A{P). As can be seen from these 
figures the Landau gauge does not give any better property compared with other 
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gauge parameters: obviously the Ward-Takahashi identity does not favor the Landau 
gauge at all ! 

§4. Conclusions and discussion 

In the present paper we solved with a numerical procedure the HTL improved 
ladder DS equation for the retarded fermion self-energy function Er to study the 
spontaneous generation of fermion mass in thermal QCD/QED, mainly focussing on 
the gauge-dependence of the solution within a general covariant gauge where the 
gauge parameter ^ is any constant number. It should be noticed that, in the DS 
equation in the point-vertex ladder approximation, no solution receives the vertex 
correction, thus the vertex renormalization constant Zi is exactly unity, Zi = 1. 
We also made an improvement in estimating the numerical integration over singular 
parts of the integration kernel in the present analysis. 

With the numerical solutions thus obtained in various gauges, we found the 
followings; 

i) In the Landau gauge the obtained solution shows a significant change compared 
to the simple ladder analyses, indicating the importance of taking the domi- 
nant effect of thermal fluctuations into the integration kernel through the HTL 
resummation procedure. 

ii) In any gauge (including the Landau gauge ^ = 0) where the gauge parameter 
^ is any constant number, the fermion wave function renormalization function 
A[P) always deviates largely from unity even at the momentum where the mass 
is defined. This fact clearly shows that the corresponding solutions explicitly 
contradict with the Ward-Takahashi identity Zi = Z2, which makes the physical 
meaning of the solution being obscure. 

iii) The obtained solutions strongly depend on the choice of gauge parameters: the 
critical temperatures (and the critical coupling constants) change significantly 
gauge by gauge. 

iv) We also determined the critical exponent v defined by Eq. (j3-ip . which controls 
how the mass i?e[M] vanishes near the critical temperature Tc- The results show 
that the temperature-dependence of mass Re[M] near the critical temperature 
Tc can be well described by the functional form Eq. (j3Tp . and that v does not 
depend significantly on the strength of the coupling, nor on the choice of gauge. 
Thus we may conclude that, despite the fact iii), the order of phase transition 
is consistent with the second order transition. 

All the above findings show the solution of the HTL improved ladder DS equa- 
tion suffers from the problem of large gauge-dependence within a general covariant 
gauge where the gauge parameter ^ is any constant number. Namely the solution 
varies significantly gauge by gauge. The most serious problem we face is there is 
no definite criterion which solution we should choose, which then reminds us of the 
fact repeatedly we have mentioned that at zero temperature any solution in the 
Landau gauge of the DS equation with the ladder kernel automatically satisfies the 
Ward-Takahashi identity Z\ = Z^'- one of the most promising criterion selecting the 
solution to have definite physical meaning. 



12 



H. Nakkagawa, H. Yokota and K. Yoshida 



Here we give some comment on the choice of the gauge in the present analysis. 
As we noted in ^ we solved the HTL improved ladder DS equation by choosing 
gauges only in the neighborhood of the Landau gauge = 0). The reason why we 
choose such gauges is as follows; i) The Landau gauge has a special significance at 
zero temperature, and might do so even at finite temperatures. With this expectation 
many analyses have been carried out in the Landau gauge, thus in performing the 
analysis in the neighborhood of the Landau gauge, we can see what really happens 
by comparing the result of our analysis with those of the preceding works, ii) There 
is a more practical reason: in our present procedure we can get nicely converged 
numerical solutions mainly in the gauges neighboring with the Landau gauge. 

Anyway with such a small change of gauge, the solution we obtained shows a 
big change in the "physical quantities" such as the critical temperature. 

Thus the only conclusion we could have from our analysis is that the chiral phase 
transition in massless thermal QED/QCD at zero fermion number density takes 
place through the second order transition. To determine the critical temperature, 
the critical coupling constant and also the corresponding critical exponents in a 
physically sensible way, we should find such a solution that satisfies at least the 
Ward-Takahashi identity Z\ = Zi- 

Needless to say, any solution of the DS equation in the ladder approximation can 
not satisfy the full Ward-Takahashi identity, stating the identity between the vertex 
function and the derivative of the fermion self-energy function. We only propose 
a possible choice of gauge where Z\ = Z2 holds at least, which may help us to 
get a physically sensible solution at the same level of significance as that at zero 
temperature analysis in the Landau gauge. It is our next plan of analysis to carry 
out this procedure and investigate the properties of such solutions. ^^^'^^^ 
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Appendix A 

Details of the numerical analysis of the DS equation 

In this Appendix, firstly we explain the problems in the numerical analysis, 
which we face in solving numerically the DS equations, Eqs. (j2-10p - (j2T2p . for the 
invariant functions A, B and C, then give the procedures we make use of to resolve 
the problems. We face two problems in the numerical analysis; 

i) As can be seen from Eq. (j2T0p . to determine the function A, firstly we perform 
the integration over the variable K in the right-hand-side (r.h.s.) of the equa- 
tion, then divide the result by p^. Needless to say, the p-dependences of both 
sides of the equation should agree with each other, and as it is easy to confirm 
that no problem appears in the analytical calculation. However, a problem does 
appear in performing the numerical calculation. In the small-p region, we are 
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forced to carry out the numerical integration in the r.h.s. in a higher accuracy 
level compared with other region of the momentum p. This is in fact a hard 
task, causing larger errors and thus being the origin of the unstable behavior 
in the numerically determined A in the small-p region. This problem becomes 
especially serious in the contribution coming from the term that depends explic- 
itly on the gauge-parameter ^, i.e., the D^^'^ term in the gauge boson propagator 
*G^'^ , Eq. ()2-2p . Formally we must divide by p^, not by p^, to determine this 
contribution to the function A. 
ii) There are several singular terms in the gauge boson propagator *G^^ , Eq. (j2-2p . 
appearing in the integration kernel of the DS equations; a) The ^-dependent 
jj^iv |;gj.jj^ pure massless double-pole mode, b) The transverse (magnetic) 
mode being proportional to B^^ receives the so-called dynamical screening, but 
it becomes massless when the energy-component of the momentum vanishes 
with the space-component of it being finite. 
In both cases we face the numerical integration over singular functions, e.g., the 
principal part and the (5- function. In the analytical calculation these functions do 
not cause any trouble, but in the numerical calculation they, especially the principal 
part, do cause troubles to obtain stable solutions. The procedures we make use of 
in order to resolve the above problems are as follows; 

i) As for the first problem i) above, we carry out the integration in the r.h.s. of 
Eq. ()2-10p by two different methods depending on the region of the momentum 
p. In the small-p region p < pth, as for the contribution coming from the D^^^ 
term that depends explicitly on the gauge-parameter ^, we expand the corre- 
sponding integrand of the r.h.s. of Eq. (|2-10p in the power series of p, keeping 
up to the p^ term, then carry out the integration. In the large-p region p > pth, 
we perform the ordinary numerical integration to get the function A. The ex- 
plicit value of the "threshold momentum" pth is determined by considering the 
stability as well as the smoothness of the solution. In the present analysis we 
choose Pth = 0.2. 

ii) As for the second problem, we use the ordinary procedure. When the integration 
over the principal part appears, we divide the integration region into two parts: 
the integration in the neighborhood of the singular point, and the integration 
away from the singular point. We can carry out the simple numerical integration 
away from the singular point. The integration in the neighborhood of the 
singular point is carried out analytically, by taking the unknown functions A, 
B, and C kept constants with the values of those at the singular point of the 
gauge boson propagator. 

In order for this method to work, the unknown functions A, B and C should 
behave smooth in the neighborhood of the singular point of the gauge boson 
propagator, and also we should take the neighborhood of the singular point as 
narrow as possible. The first point can be checked a posteriori by the obtained 
solutions, and the result is satisfactory. As for the size of the neighborhood 
where the integration is performed analytically, we must choose by seeing the 
stability and the smoothness of the obtained solutions. In the present analysis, 
we choose \po — /cqI < 0.04 and |p — k| < 0.2. 
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